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Abstract 



^ We derive the extension of the Hubbard model to include Rashba spin-orbit coupling that correctly describes Aharonov-Bohm and 



Aharonov-Casher phases in a ring under applied magnetic and electric fields. When the ring is connected to conducting leads, we 
develop a formalism that is able to describe both, Kondo and interference effects. We find that in the Kondo regime, the spin-orbit 
H coupling reduces strongly the conductance from the unitary limit. This effect in combination with the magnetic flux, can be used 
io produce spin polarized carriers. 
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1. Introduction 

There is a great interest in nanoscopic systems, either be- 
cause of its potential application in nanodevices or as ideal 
systems to test theories for non-trivial physical problems. 
In particular, the Kondo effect is present in many of these 
systems, like magnetic impurities on clean noble metal sur- 
faces or quantum corrals p] , small clusters on these surfaces 
|2|3j and systems of quantum dots |4|5|6|7|8(9] . In particu- 
lar, the unitary limit, which consists in the maximum pos- 
sible conductance through a quantum dot has been reached 
experimentally [4]. The Kondo model out of equilibrium 
brings new challenges to the theory [10|llj . 

On the other hand, effects of interference in quantum 
paths and the Aharonov-Bohm effect have been demon- 
strated in mesoscopic rings with embedded quantum dots 
|4I5I12I13| . 

The calculation of transport through a mesoscopic ring 
in which both interference effects and interactions leading 
to the Kondo effect are present is not trivial. An example 
of this is a recent debate about the role of interactions in 
dephasing [14] . Even knowing the exact eigenstates of the 
ring, there is no simple procedure to calculate the conduc- 
tance. When the coupling V of the ring to the conducting 
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leads is small, using perturbation theory up to second order 
in V , an expression has been derived, which is also exact for 
any V in the non-interacting limit, the Jagla-Balseiro for- 
mula [15] . Similar equations were used recently |16I17| . Al- 
ternative perturbative expressions were also proposed [18] . 
Unfortunately these expressions are not valid in the Kondo 
regime, in which the ground state of the isolated ring is 
Kramers degenerate (odd number of electrons), because 
the ground state of the whole system is a singlet with a 
characteristic energy scale (the Kondo temperature) Tx ~ 
W exp[—l/p J], where W is the band width, p the den- 
sity of states at the Fermi level and J rj V 2 [see Eq. (TT2"]) or 
Ref. 3]], which cannot be recovered by perturbation the- 
ory in V. Previous calculations of the conductance through 
strongly correlated rings in which the effects of interference 
were important to detect signatures of spin-charge sepa- 
ration, assumed that a Zeeman term destroys the Kondo 
effect in the system [16]. For a ring described by the ionic 
Hubbard model, it has been shown that the conductance 
through the system is related to the quasiparticle weight 
[17] . The physics of the Kondo regime was recovered by 
mapping the model into an impurity Anderson model, but 
at the cost of losing interference effects. 

One of the branches of nanophysics with particular re- 
cent interest is spintronics, which consists in developing 
means of creating and controlling spin polarized currents 
in nanoscale systems [19] . An important ingredient for this 



purposes is the Rashba spin-orbit coupling 

H S oc = aa-E x (p-eA), (1) 

which is present in quantum wells. Recent experiments in 
semiconductor mesoscopic rings have shown that the con- 
ductance oscillates not only as a function of the applied 
magnetic field (Aharonov-Bohm effect) but also as a func- 
tion of the applied electric field E perpendicular to the plane 
of the ring (Aharonov-Casher effect) [21122) . In this effect, 
the electrons, as they move, capture a phase that depends 
on the spin, as a consequence of spin-orbit coupling. The 
main features of the experiment can be understood in a 
one-electron picture 23 24J. However, this picture is inad- 
equate in the presence of strong correlations. In particular, 
the Kondo effect cannot be described. 

The aim of the present work is two-fold: first, to pro- 
vide a method to calculate the conductance through a ring 
that is able to capture at the same time the Kondo physics 
and interference effects; second, to apply the formalism to 
calculate the conductance through a ring under the action 
of both, magnetic and electric fields together with Rashba 
spin-orbit coupling i?soc- I n the second task we stum- 
bled with the difficulty of finding the form of the Hub- 
bard model in the presence of Hsoc, which adequately de- 
scribes the periodicity with the applied electric field seen 
in the experiments. We solved this problem using a non- 
abelian gauge transformation in the continuum version of 
the model, which translates the effect of Hgoc to a change 
in the boundary conditions. 

The main physical result is that the combination of both 
Kondo effect and spin-orbit coupling leads to a strongly 
spin dependent conductivity that might be used for spin fil- 
tering purposes. Some results have been recently published 
[25] . In section 2 we briefly explain the general method to 
calculate the conductance. In section 3, the Hubbard model 
including the effect of the spin-orbit coupling Eq. |T]) is de- 
rived. Section 4 contains results for a four-site ring. Section 
5 is a summary and discussion. 

2. Formalism to calculate the conductance 

To calculate the conductance through an interacting sub- 
system, it is necessary to obtain some Green functions of 
the whole system, which includes the interacting part both 
leads 26J . An example of an interferometer assembled ex- 
perimentally with one quantum dot is given in Fig. 1 of 
Ref. 5. The Hamiltonian can be written as 

H = H I + H l + H v 
Hi = B'v -VgJ^dUi*, 

icr 

— QO OO 

Hi = t c { 4-i^ c i<? + 4+i,* c i<? + H - c )> 

t — 0,cT i — 1,CT 

H v = V(J2 c Ldoa + c\j la + H.c). (2) 



Hj describes the interacting subsystem (in the case de- 
scribed below, a Hubbard ring with spin-orbit coupling). 
The operator d ia creates an electron with spin a at site i 
in this subsystem. The second term of Hi represents the 
effect of a gate voltage. The term Hi describes both leads, 
and Hy is the coupling between both leads and the inter- 
acting region. We label the sites in this region so that site 
is coupled to the left lead and site 1 to the right one. The 
prime in H{j reminds us that the Hamiltonian of the sub- 
system is written in a particular gauge, as described in the 
next section. Note that the phases of all Ci„ operators can 
be chosen so that t c and V are real, regardless of the above 
mentioned gauge. The Zeeman term is neglected here. To 
reduce the number of parameters, we have assumed identi- 
cal left and right leads and the same coupling between any 
of them and the interacting subsystem. 

Since Hi describes a finite system, it can be diagonalized 
exactly. Our general approach to calculate the conductance 
amounts to a truncation of the Hilbert space of Hi , retain- 
ing only two neighboring charge configurations, with n and 
n — 1 particles. It can be shown that this procedure is valid 
for small enough V [5] • Calculating the matrix elements of 
Hy in the truncated Hilbert space leads to a generalized 
Anderson model 

ff G AM=if ; +X>i^")(V>"| 

+v ^wiAM-'mi+K-c.), (3) 

where [ip?) and E™ denote the j-th eigenvector and eigen- 
value of Hi in the configuration with n particles and 

Pl 3 « = (r k - l \d m \^) (V = 0,1). (4) 

In the simplest case, only the spin singlet (doublet) 
ground state of Hi is relevant for the configuration with 
an even (odd) number of electrons and Hgam reduces 
to the ordinary Anderson model [17]. However, interfer- 
ence effects, for example the vanishing of the conductance 
through the ring for certain values of the flux [16| is a con- 
sequence of an orbital degeneracy of levels for that flux and 
more states than one Kramers doublet should be included 
in Hqam, as we report in section 4. 

In general, Hqam can be represented using one or more 
slave bosons that represent the relevant eigenstates of Hi 
with an even number of particles [27)28] . A model with a 
doublet hybridized with a singlet and a triplet has been 
solved using the numerical renormalization group [29]. 
From the solution of Hq am, the conductance is obtained 
using known expressions that relate it with the exact Green 
functions of Hqam [2S] • 

3. Hubbard model in the presence of spin-orbit 
interaction 

The simplest way to add the physics of the spin-orbit cou- 
pling in the Hubbard model seems to be to replace -ffsoc 



given by Eq. (|TJ) by a tight-binding version in which the 
current in each link is multiplied by the Pauli matrix per- 
pendicular to the link and to the electric field, to build the 
cross product of Eq. ([T]) [30] . For an electric field in the z 
direction and a ring in the x, y plane this gives: 

H tb = aE z h/(2a) ^[icos <p(d' i+n dii + 

i 

+ sin (p(d\ +n d iL - dl +n dii) +H.c.]. (5) 

However, in the experimentally assembled rings it is clear 
that the conductance oscillates as E increases [21122] and 
therefore one expects that the electric field enters an expo- 
nential as a phase, and this is not apparent in Eq. (|5|). 

To follow a procedure that leads to such an exponential 
dependence, we consider the version of the model in the 
continuum Ufj- This is obtained from the noninteracting 
version derived by Meijer et al. [31j for electric and magnetic 
fields perpendicular to the plane of the ring, adding a local 
interaction. The Hamiltonian is 
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m 



. d 7 

~ l a X~ + o^Wj, 

dtpj 0o 2 



(6) 



l<3 



where 



is the azimuth of the 7-th electron, Q = 



h/(2m*r ), to* is the effective electron mass, r is the ra- 
dius of the ring, 7 = aE z /{rQ.) is proportional to the 
Rashba constant a and the electric field E z , <fr — Birr 2 
is the magnetic flux, 0o = hc/e is the flux quantum and 
<j r (ip) — a x cos ip + a y simp is the Pauli matrix in the 
radial direction. 

It can be easily checked that the unitary transformation 
T = IL*(<Pj)> with 



t(ip) = exp 
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n$ = (— sin 6, 0, cos 6), 6 = arctan (7), 



ip' = ip\/l + ~f 2 , 



(7) 



makes the dependence of the fields [0 and 7 in Eq. (J6j> ] 
disappear in the transformed Hamiltonian — T^HfjT . 
Therefore H£} is the continuum version of the ordinary Hub- 
bard model, in which the magnetic flux and the spin-orbit 
coupling have been gauged away. The price to pay is that 
the transformed one-particle spinors satisfies the bound- 
ary conditions x'{2tt) = v(2tt)x'{Q) (instead of periodic 
ones). This implies that in the Hubbard model for a ring 
of N sites, the last hopping (from angle (p = —2ir/N to 0) 
should be modified. Diagonalizing T^(2ir) one obtains the 
eigenvalues exp[i(<&AB + o~&ac)]i where a = ±1 for spin 
pointing in the ±ng direction, &ab — ^4 > /4'o, and $ac = 

7r{[l+7 2 ] 1/2 -l}- 

Then, the transformed Hubbard model for a ring of TV 
sites takes the form 



H 'u = - E t[4+i,di,+H. 

i=0,a 

-t [e^^+^^d^djv-i,* + H.c. 



where 



4> 



AC 



= Vir 2 + R 2 -tt, 



(8) 



(9) 



and the ratio R — tt^ is given by HaE z N/ (2ta) if the mass 
is eliminated from the curvature of the bottom of the non- 
interacting band (m* = ft 2 /(2ta 2 ), where a — 2irr/N is the 
lattice parameter). If instead to* — hkp/vp is used, where 
kp, Vp are the Fermi wave vector and velocity respectively 



R 



HaE z Nkpa 
2ta sin(fc^a) 



(10) 



In any case, the same qualitative physics is described. The 
rotational invariance of H, 



u can be recovered replacing 
d] a = d] a exp[i(<S> AB + a<f> AC )j/N}. 

It is important to remark that the spin quantization axis 
of Eq. ((5|) varies with position. For <p = (corresponding to 
y = 0), it is given by fig (see Eq. ([7])) and it lies in the x, z 
plane. For other angles the spin quantization axis rotates 
with ip in such a way that the component in the plane of 
the ring (cc, y) always points towards the center of the ring. 

It is interesting to note that Berry phases captured in 
adiabatic evolutions of a system in &ab and &ac give in- 
formation on the polarization, opening of a spin gap, phase 
transitions and ferrotoroidic moments [32 33J . 

4. Results for 4-site ring 

For the explicit calculation, we take a ring of N = 4 sites, 
connected to the leads at opposite sites, lying at (f> = and 
cj) — 7r. We assume that the leads are described by a con- 
stant density of states po = 1/W and set the band width 
W = 60£ (much larger than the hopping t in the ring) . The 
Fermi energy of the leads is set at the on-site energy in the 
ring (zero) .We have included all doublet states with n = 3 
and the singlet ground state for n = 4 of Hj, to construct 
the generalized Anderson model if gam (see Eq. ©). We 
have solved Hqam in a slave boson mean-field approxima- 
tion, which is known to reproduce correctly the exponential 
dependence of the Kondo scale on the parameters. Details 
are given in Ref . [25] . 

4.1. The different regimes 

The properties of the effective generalized Anderson 
model, Eq. (J3j) and therefore the conductance through the 
ring, differ according to different regimes which depend 
on the ratio of the charge-transfer energy E^ — E^\ 
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Fig. 1. Different regimes of the model for spin-orbit strength R = 
and two values of the applied magnetic flux:$Afl = (top) and 
<J>AB = 1" (bottom). 

f n) 

where Eq is the ground state of Hi in the subspace of n 
electrons, and the effective resonant level width 

A = n P0 V 2 J2\^\ 2 - (11) 



If the four-particle singlet is well below the lowest 
Kramers doublet (E^ - E^ ] > A) the syste m is m 
the non- magnetic regime. For increasing values of E^ , 
the system enters first the intermediate valence zone 



\E 



(3) 



E, 



(4) 



A), and then the Kondo regime when the 
lowest spin doublet is well below the sing let (E^-E^» 
A). 

In Fig. Q] we show the different regions of parameters of 
the model U, V g and V (setting t as the unit of energy) cor- 
responding to the different regimes in absence of spin-orbit 
coupling. The boundary between the mixed-valence regime 
and the nonmagnetic (Kondo) one has been defined arbi- 



trarily in the figure by E { 



(3) 
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E, 



(4) 



-8A). 



with the applied magnetic flux. The four-particle singlet is 
favored for &ab = tt. This can be understood already for 
the non-interacting ring, since the ground-state energy for 
four particles E^ is optimized for <&ab = 7T, while it passes 
through a relative maximum at &ab = 0. In the strongly 
interacting case U — * +oo, Hjj reduces to a t — J model 



and Eq 4 ^ becomes independent of the flux, while Eq 3 ' is 
minimized for $ab =0 16J. 

Clearly, there is no mixed valence regime for V = 0, 
since this implies A = 0, while the extension of this regime 
increases nearly quadratically with V. 



(3) 



4.2. The non-magnetic regime 
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2. Conductance as a function of magnetic flux for R = 0, 
— 0.8t, U = 2t and several values of V. Full lines: our formalism. 
Dashed lines: Jagla-Balseiro formula. 

In Fig. [2] we show the conductance G — Gf + as a 
function of magnetic flux in the non-magnetic regime for 
different values of V and without applied electric field (R = 
0). The quantum of conductance is denoted as Go = 2e 2 /h 
. The total conductance is even with flux (f> and we show G 
in the interval < (j) < 4> /2 (or < $ab < 7r). 

The conductance is very small due to the fact that there 
are no available states of the ring near the Fermi energy. G 
increases with the hopping to the leads V, since this term 
promotes electrons to the leads at the Fermi level. In addi- 
tion, in this regime correlations play a minor role and one 
expects that the Jagla-Balseiro formula [TS], which is exact 
in the non-interacting case, gives accurate values for the 
conductance. In fact, our results show the same qualitative 
behavior and for small V it can be demonstrated that both 



i 



approaches are equivalent in this regime. However, there 
are significant quantitative differences for V/t = 1. It is dif- 
ficult to state which approach is the most acurate in this 
case. 

In absence of spin-orbit coupling, for an applied flux of 
half a flux quantum, the conductance vanishes due to de- 
structive interference, as a consequence of reflection sym- 
metry of the ring for a non-degenerate ground state [16] . 

4.3. Spin dependent conductance in the Kondo regime 

In the rest of this paper, we consider Eq 3 J < , when 
the ring is the mixed valence or Kondo regime. In Fig. [3] 
we show the conductance as a function of the applied mag- 
netic flux fyis for several values of the applied electric field 
E z . We discuss first the case E z = R = 0. In contrast to 
the results in the non-magnetic regime, the conductance 
takes appreciable values, being near to the ideal one for 
0.2 < &ab < 0.5 and 1.5 < $ab < 1-8. The system can 
be considered to be in the Kondo regime for these values 
of Q ab and in the mixed valence regime for other fluxes. 
This is also consistent with the calculated occupation of 
the levels, shown elsewhere |25] . which is nearly one in the 
Kondo regime. As expected from Fig. [2 the region of $ab 
for which the system is in the Kondo regime increases with 
decreasing V, and for V — 3.1t nearly perfect conductance 
is obtained for 0.1 < $ab < 0.7 and 1.3 < $ab < 1.9 [25] . 
This corresponds to the unitary limit. This limit, observed 
experimentally in systems with one quantum dot different 
than ours [4], is characteristic of the Kondo regime. The 
Jagla-Balseiro formula fails to describe this case, giving 
very small values of the conductance. Note that the con- 
ductance still vanishes at <&ab = 7r as a consequence of 
interference. This together with the ideal conductance for 
the Kondo regime shows that our approach includes both, 
interference and Kondo effects, in contrast to previous cal- 
culations |16I17| . 



Precisely in the Kondo regime, where the conductance 
reaches its maximum possible value, is where the effect of 
the spin-orbit coupling times the applied electric field E z 
is most noticeable. As it can be seen in Fig. [31 moderate 
values of E z are enough to decrease the conductance from 
Go to very small values. The main reason for this is that 
E z affects both spin projections in opposite ways breaking 
SU(2) symmetry, and therefore tends to destroy the Kondo 
singlet. The stabilization energy of this singlet is of the 
order of the Kondo temperature 



T K ~ W exp{-Tr(E { 4) 



4 3) )/A), 



(12) 




and therefore decreases exponentially inside the Kondo 
regime. For the parameters of Fig. [31 the system is deeper 
in the Kondo regime for $ab ~ ±0.37r (mod 2tt). Due to 
the loss of reflection symmetry, in presence of spin-orbit 



Fig. 3. Conductance as a function of flux for V = 3.5t,V g = 0, U = 6t 
and different values of the electric field times spin-orbit coupling 
strength R (see text). 

coupling, complete conductance cancellation is not realized 
at = tt. 

The above mentioned breaking of spin SU(2) symmetry 
as a consequence of the spin-orbit coupling leads to signifi- 
cant differences in the conductance G a for both spin orien- 
tations a. Since in our Hamiltonian, the quantization axis 
depends on the position of the electron, this conductance 
should be interpreted in the following way: if a spin a (up or 
down) in the quantization direction ng = (— sin#, 0, cos#) 
(see Eq.([7])) is injected in the ring at the right lead (ip = 0) 
it comes out at the left lead (tp — tt) with spin a in the di- 
rection fig = (sin#, 0, cos 6) or vice versa. For an arbitrary 
direction, the incident wave should be decomposed in the 
components along the corresponding quantization axis. In 
Fig. [4] we show the degree of polarization of the conduc- 
tance P = (Gj — GjJ/(Gj + GjJ as a function of &ab for 
several E z . The ratio P = ±1 at symmetric points near 
^ab — Ti" because one of the G a vanishes there and the 
other one at the symmetric point around n [the Hamilto- 
nian is invariant under time reversal and change of sign 
of $ab = 2(7rr) 2 S/(/)o]. For fluxes corresponding to the 
Kondo regime P ~ 0.4 and one of the G a can still be near 
the ideal one. 



5. Summary and discussion 

We propose an approach to calculate the conductance 
through a ring of interacting quantum dots, weakly cou- 
pled to conducting leads, that takes into account both in- 
terference effects and non-perturbative many-body ones, 
mapping the relevant states into an effective generalized 
Anderson model containing several non-degenerate states. 
Here we have solved the model using a slave-boson repre- 
sentation in the saddle-point approximation, but alterna- 
tive treatments like the numerical renormalization group 
are also possible. 

We have derived an extension of the Hubbard model 
that includes spin-orbit interaction, absorbing it in op- 
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posite Aharonov-Casher phases for spin up and down in 
an adequately chosen quantization axis, defined by a non- 
abelian gauge transformation. For a non-interacting sys- 
tem, it has been noticed previously that the Rashba spin- 
orbit coupling can be gauged away in a ring 34J , although 
the explicit dependence of the boundary conditions [the 
Aharonov-Casher phase $>ac, see Eq. ((9])] with the applied 
electric field was not given. Using Bethe ansatz, the persis- 
tent currents in a Hubbard ring were calculated as a func- 
tion of &ac by Fujimoto and Kawakami [35]. These authors 
also note that the trick used to gauge away the spin-orbit 
coupling can be extended to any SU(2) invariant interac- 
tions (not only local as U). For example a nearest-neighbor 

repulsion Yji aa > d\+ia d i+^A^^<y' is not modified by a spin 
rotation at any site. Therefore, the thermodynamic prop- 
erties of the extended Hubbard model including this inter- 
action should not vary with spin-orbit coupling in the ther- 
modynamic limit, in which boundary conditions are irrele- 
vant. This fact is not obvious in alternative treatments [36] . 

We have calculated the conductance through a ring de- 
scribed by the Hubbard model with Rashba spin-orbit cou- 
pling in presence of magnetic and electric fields. The main 
effect of the spin-orbit coupling is to tend to destroy the 
Kondo effect, leading to a strong spin dependence of the 
conductance. 
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